Coherent instabilities in a semiconductor laser with fast gain recovery 
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We report the observation of a coherent multimode instability in quantum cascade lasers (QCLs), 
which is driven by the same fundamental mechanism of Rabi oscillations as the elusive Risken- 
Nummedal-Graham-Haken (RNGH) instability predicted 40 years ago for ring lasers. The thresh- 
old of the observed instability is significantly lower than in the original RNGH instability, which we 
attribute to saturable-absorption nonlinearity in the laser. Coherent effects, which cannot be repro- 
duced by standard laser rate equations, can play therefore a key role in the multimode dynamics of 
QCLs, and in lasers with fast gain recovery in general. 

PACS numbers: 42.55.Px, 42.60.Mi, 42.65.Sf 



The fundamental coherent mechanism that can desta- 
bilize a single-mode laser was predicted in the early 60s 
l] and was later extended to multi-mode lasers [2, Q 
where it became known as the RNGH instability. These 
instabilities became classic landmarks for the general 
field of nonlinear dynamics 0, Q because they emerge 
in conceptually the simplest laser model, which in the 
single-mode case was shown to be equivalent to the 
Lorentz model of deterministic chaos Q . Another feature 
that makes these instabilities so interesting and unique is 
their coherent nature that involves the polarization of the 
medium as a nontrivial dynamical variable. Most other 
physical mechanisms that can drive a laser from a single- 
mode to a multi-mode regime, such as spatial and spec- 
tral hole burning, Q-switching, and saturable absorption 
0, S] , can be adequately described within the standard 
rate equation formalism, in which the polarization of the 
active medium is adiabatically eliminated. The RNGH 
instability and its single-mode twin cannot be explained 
by the rate equations. Such coherent effects can be only 
observed when the polarization is driven faster than or 
comparable to the dephasing time T2 Q. 

The origin of both phenomena is the oscillation of the 
population inversion at the Rabi frequency f^Rabi that 
takes place when the intracavity laser intensity becomes 
large. This results in a modification of the gain spectrum 
and the emergence of sidebands separated from the max- 
imum of the gain curve by an amount corresponding to 
the Rabi frequency. These sidebands can be regarded as 
a manifestation of parametric gain. The instability sets 
in when the intracavity power is sufficiently large: the 
Rabi angular frequency f2R a bi has to be greater than the 
relaxation timescales of the gain medium (more precisely, 
^Rabi is sufficiently greater than (TiT 2 ) _1 / 2 , where T\ is 
the gain relaxation time). The instability threshold is of- 



ten called the second laser threshold due to its universal 
nature. 

Pioneering theoretical works stimulated extensive ex- 
perimental studies that finally resulted in the observa- 
tion of the Lorenz-type chaos in a far-infrared single- 
mode laser [lfj]. However, despite almost 40 years of ef- 
forts, the experimental demonstration of the multimode 
RNGH instability has remained somewhat controversial 
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In lasers with long gain recovery compared to the cav- 
ity roundtrip time, the instability caused by a saturable 
absorber can often lead to mode locking 8]. When the 
gain recovery time is short compared with the cavity 
round-trip time, it is usually assumed that laser dynamics 
becomes very primitive and uninteresting (so-called class 
A laser). In this case mode locking is impossible accord- 
ing to conventional theory, and the relaxation oscillation 
frequency becomes purely imaginary }17l |. Surprisingly, 
as we show in this Letter, it is under these conditions 
that the elusive RNGH instability can be observed. We 
show that quantum cascade lasers (QCLs) are uniquely 
suited for studying these coherent effects which, along 
with spatial hole burning (SHB), become a key factor in 
dictating the dynamics of the laser. 

QCLs, because they are based on ultrafast tunneling 
and phonon-limited intersubband transitions, belong to 
the class of lasers which have a extremely fast gain recov- 
ery, in the range of a few picoseconds [l8( ■ Recent exper- 
iments showed indeed that the gain recovers within a few 
picoseconds, which is approximately an order of magni- 
tude shorter than the cavity round-trip time [lj|. Since 
its invention in 1994, QCLs have undergone tremendous 
improvement 12 CM . Recent development of low loss, high 
power QCLs [2lL I22] ] enables the study of previously 
under-investigated aspects, such as the richness of the 
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FIG. 1: (a) Optical spectra vs. pumping ratio (j/jth) above threshold obtained in cw at 300K with a 3 /im wide buried 
heterostructure lasers emitting at 8.38 /im. For 1 < j/jth < 1-2 the spectra are identical to j/jth < 1-2. (b) Spectral splitting 
and twice the Rabi frequency fiR a bi/(27r) vs. square root of output power collected from a single laser facet. The different 
quantities reported on the graph were deduced from the experimental data shown in [T] (a) . The dashed line is a least-square 
linear fit of the data. 



optical spectrum and the ultrafast dynamics of these de- 
vices. In Ref. 23], strong evidence of self-pulsations at 
the cavity roundtrip frequency was reported in QCLs, in 
particular a large broadening of the spectrum above the 
threshold of this instability was observed. However, no 
detailed pulse characterization was provided. The tech- 
nological potential of QCLs calls for a better understand- 
ing of the interplay of various instabilities in the param- 
eter regime dictated by these lasers. Moreover, the Rabi 
frequency in QCLs at the power levels of a few hundred 
milliwatts is of the order of a few THz, much larger than 
the spacing of Fabry-Perot modes. Therefore coherent 
effects should be easily observable in QCLs. 

In this Letter we present a clear experimental demon- 
stration of a coherent instability, driven by the same 
mechanism as the RNGH instability. It is identified in 
the most direct manner, by demonstrating in the optical 
spectrum of QCLs a splitting corresponding to twice the 
Rabi frequency. To the best of our knowledge, this is 
the first observation of the RNGH mechanism in a semi- 
conductor laser and even more generally in a solid-state 
laser. 

The instability observed differs in some respects from 
the original RNGH instability 0, 0|. The threshold of 
instability is as low as ~ 50% above the laser thresh- 
old. In addition, the pure RNGH instability typically 
gives rise to spectra with one central mode and two side- 
bands separated from it by the Rabi frequency, whereas 
in our experiments we observed two peaks only, similarly 
to Ref. ll| . However the mechanism of the instability is 
the same in essence, namely the Rabi oscillations of the 
population inversion due to coherent effects. The dif- 
ferences from the RNGH instability as it occurs in ideal 



conditions 0, 01 can be attributed to the presence of sat- 
urable absorption and SHB. 

The QCLs studied were fabricated from two different 
wafers (wafer # 3251 and 3252) grown by metalorganic 
vapor phase epitaxy. The devices were processed into 
buried heterostructure lasers, in which an insulating Fe- 
doped InP layer is regrown after etching of the ridges 
[2l[ |22 |. The active region of all the samples tested is 
based on a four-quantum-well design, which rely on a 
double phonon resonance to achieve population inver- 
sion [24J. Note however that the multimode operation 
described in the present letter was also observed with 
lasers based on so-called three-quantum- well designs [l8| . 
Fig. QTa) shows the optical spectra of a laser operated in 
continuous wave (cw) at room temperature. The active 
region of this laser is 3/xm wide and its emission wave- 
length is close to 8.38/^m (wafer #3251). The laser was 
cleaved into a 3.25mm long bar and soldered with Indium 
onto a copper heat sink. The spectra were measured 
by a Fourier transform infrared spectrometer (FTIR) 
equipped with a deuterated triglycine sulphate (DTGS) 
detector. 

As shown in Fig. [T] (a), the laser spectrum is single 
mode close to threshold and broadens as the pumping 
current increases, splitting into two separated humps. 
The difference between the weighted centers of the two 
peaks increases linearly as a function of the square root 
of the collected output power from one facet, as shown 
in Fig. [T] (b) (square dots with the dashed line as its 
best- fit). The Rabi angular frequency f^Rabi can be 
easily calculated using the formula f^Rabi = [lE/H = 
fiy / 2nl ave /(ce)/h, where \i is the electron charge times 
the matrix element of the laser transition (=2.54nm). 
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iavo is the average intracavity intensity in the gain re- 
gion, which can be derived from the measured output 
power, c is the speed of light in vacuum and n is the 
background refractive index. For all the values of the in- 
tensity corresponding to the spectra reported in Fig.QTa), 
^Rabi was calculated, multiplied by a factor two and then 
added to Fig.[2](b) (solid line). A very good agreement is 
found between the experimental splitting and twice the 
estimated Rabi frequency. Both curves fall indeed well 
within the error bars j25|. As mentioned before, the the- 
ory behind the RNGH instability predicts that the large 
intracavity intensity will result in parametric gain at fre- 
quencies detuned from the maximum of the gain curve 
by the Rabi frequency. The agreement mentioned above 
is thus a strong indication of the RNGH instability in 
QCLs. 

In order to better understand the experimental spec- 
tra of the QCLs presented in Fig. [T] (a), we use a simple 
model based on the standard one-dimensional Maxwell- 
Bloch equations Q , where the active medium is described 
by an "open" two level system [26| . However contrary to 
the standard unidirectional Maxwell-Bloch equations, we 
allow the electromagnetic field to propagate in both di- 
rections. The waves traveling in the two directions are 
coupled, as they share the same gain medium. This gives 
rise to SHB [7(: The standing wave formed by a cavity 
mode imprints a grating in the gain medium through gain 
saturation. As a result, other modes become more favor- 
able for lasing, and a multimode operation is triggered. 

In the slowly varying envelope approximation, the 
equations read: 
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The + and — subscripts label the two directions of prop- 
agation. E and r) are the slowly- varying envelopes of the 
electric field and the polarization respectively. The actual 
electric field and polarization are obtained by multiplying 
E and ry by e lut (u> is the optical resonance frequency) and 
taking the real part. The position-dependent inversion 
is written as the sum of the three terms, Ao, A^e 2lfez , 
and A^~e _2jfcz , where (A^~)* = A^. The inversion is 
thereby represented by two slowly varying functions (Ao 
and A%), and e ±2lkz gives the fast variation in space. All 
the quantities mentioned so far are functions of space z 
and time t. 

£(E + ,E^) is the loss in the cavity (not including the 
mirror loss) , which is allowed to be nonlinear and to de- 
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FIG. 2: Results of numerical simulations of the spectra based 
on the Maxwell-Bloch equations including a saturable ab- 
sorber and spatial hole burning for different values of the 
current density normalized to the threshold value. 



pendent on the intensity. In this work we assume 

£(E + ,E_) = l - 1 (\E + \* + \E^), (2) 

where £q is the linear loss and 7 is the self-amplitude 
modulation coefficient characterizing the nonlinear (sat- 
urable) part of the loss. Such a saturable absorbtion 
mechanism can come from Kerr-lensing 23 1, caused 
by a nonlinear refractive index rii in the active region. 
As the intensity increases, the mode is more confined in 
the plane transverse to the propagation direction, and the 
net gain it experiences is greater. The reason is twofold: 
First, the mode overlaps more with the active region, 
leading to a larger modal gain (this mechanism is often 
called "soft Kerr-lensing" [27|). Second, the overlap with 
the metal contacts is smaller, leading to smaller losses 

E+ and £7_ satisfy the boundary conditions E+ = rE- 
at the z — boundary and rE + = E- at the z = L 
boundary (L is the cavity length and r « 0.53 is the re- 
flection coefficient). The other quantities in Eq. ([I]) are 
constants: k, N, and T are the wavenumber (in the ma- 
terial) associated with the resonance optical frequency, 
the electron density in the active region, and the over- 
lap factor between the laser mode and the active region 
respectively. 

Figure [2] shows spectra that were obtained by solv- 
ing numerically the equations Eq. |T]) with the follow- 
ing parameters: for the saturable absorber, we used 



7 : 



10 8 ^§, obtained from two-dimensional mode simu- 



lations, assuming a ni ~ 10~ 9£ ^- [23j. The index change 
due to this 712 at typical intracavity intensities is about 
10 -3 . The other parameters are £q = 5cm -1 , T\ = 0.5 
ps [lj|, T 2 = 0.07 ps (corresponding to a gain FWHM 
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bandwidth of 4.8 THz), L = 0.3 cm, and n = 3, which 
are typical values for these lasers. N and T are not needed 
as long as the pumping is expressed relative to the las- 
ing threshold. Note that the simulated spectra presented 
in Fig. [2] are averaged over about a microsecond. Only 
then does the average spectrum reach a steady state and 
a clear pattern shows up. The averaging is motivated 
by the fact that experimentally the spectra are acquired 
over an even much longer timescale. The envelopes of 
the spectra show two clear peaks whose separation com- 
pares well with twice the Rabi frequency, similarly to the 
experiment. 
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FIG. 3: Optical spectra vs. pumping ratio above threshold 
obtained in cw at 300K with a 5 fim wide buried heterostruc- 
ture lasers emitting at 5 /jm. 

The lowering of the RNGH instability threshold by 
a saturable absorber can be established analytically by 
means of linear stability analysis. We propose this mech- 
anism as the main reason for the observation of the 
RNGH instability at much lower pumping than RNGH 
theory predicts. In order to support this idea, we 
now present spectra from another device similar to the 
one described previously. The only difference between 
the two lasers is a shorter optical wavelength (5/im) 
(wafer #3252) and a wider active region (5/um). The 
two-dimensional waveguide simulations indicate a much 
weaker Kerr-lensing effect in these QCLs (7 is smaller by 
a factor of 4) , due to the much larger ratio of active re- 
gion width to wavelength. The measured optical spectra 
obtained at 300K in cw mode with the 5/im device are 
shown in Fig. [3J The data clearly show that the laser is 
at first single mode close to threshold and becomes multi- 
mode immediately after a slight increase of the pumping 
current. The envelopes of the spectra consist of multiple 
peaks, with an average separation 0.2THz, independently 
of the pumping. Numerical integration of Eq. (JTJ with- 
out a saturable absorber (7 = 0) leads to spectra that 



qualitatively agree with the ones in Fig. [3] 

Ref. [15] suggested that the suppression of the central 
peak in RNGH-type spectra can be due to the complex 
level structure of the gain medium, a dye molecule in that 
case. We show that SHB can also result in the suppres- 
sion of the central peak (Fig. 

Our postulation of saturable absorption due to Kerr- 
lensing is supported by more extensive study of different 
devices beyond those shown in this Letter. First, we ob- 
served that for the same emission wavelength, a broad ac- 
tive region leads to a less pronounced RNGH-type signa- 
ture. Second, we have also tested several standard ridge 
waveguide QCLs, for which the sidewalls of the ridges are 
covered by the gold contact. For these devices the cou- 
pling between the optical mode and the metal is expected 
to be stronger and so is the effect of saturable absorber 
due to Kerr-lensing. The spectral behavior observed in 
this class of devices is dominated by RNGH-type insta- 
bility. 

In summary, a coherent multimode instability in quan- 
tum cascade lasers (QCLs) has been observed. It is 
similar in many ways to the Risken-Nummedal-Graham- 
Haken (RNGH) instability. The threshold of the ob- 
served phenomenon is significantly lower than in the orig- 
inal RNGH instability, which is attributed to the pres- 
ence of a saturable absorber in the laser. For devices 
with a weaker saturable absorber, the envelope of the 
optical spectrum consists of many maxima whose sepa- 
rations are independent of the intracavity power. The 
nontrivial shape of the spectrum can be explained by 
SHB. 
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